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The Lieb-Robinson theorem states that the speed at which the correlations between two distant 
nodes in a spin network can be built through local interactions has an upper bound, which is called 
the Lieb-Robinson velocity. Our central aim is to demonstrate how to observe the generalized 
Lieb-Robinson velocity in an Ising spin chain with a strong transverse field. Here the generalized 
Lieb-Robinson velocity is defined as the correlation propagation speed from a given initial non- 
correlated state, which is upper bounded by the Lieb-Robinson velocity. We adopt and compare 
four correlation measures for characterizing different types of correlations, which include correlation 
function, mutual information, quantum discord, and entanglement of formation. All the information- 
theoretical correlation measures demonstrate the existence of the generalized Lieb-Robinson velocity. 
In particular, we find that there is a sudden switch of the generalized Lieb-Robinson speed with the 
increasing of the number of spin. 

PACS numbers: 03.65.Ud, 03.65. Vf, 03.67.Mn 
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I. INTRODUCTION 

The correlations in a quantum many-body system 
make its physical properties can not be regarded as the 
simple sum of the properties of the composite subsys- 
tems. Recent studies of correlations in quantum infor- 
mation science show that correlations can be used as in- 
dispensable resources in completing some computing and 
information tasks Thus the correlations play a key 
role both in quantum many-body physics and in quantum 
information science. 

The correlations in a many-body system can be classi- 
fied into different types according to different standards, 
where the two-party correlation is relatively well under- 
stood and widely used in practical problems. The two- 
party correlations in a many-body quantum state are 
measured by correlation functions associated with phys- 
ical observables in traditional physics, while they are 
measured by the mutual entropy in quantum informa- 
tion. The two-party correlations can be further classified 
into quantum two-party correlations and classical two- 
party correlations, and entanglement is a specific type of 
quantum correlations, which is extensively investigated 
in quantum information science. The measures of differ- 
ent types of corrrelations are proposed in literature. The 
widely used entanglement measure is the entanglement of 
formation (EoF) 0,3, and quantum discord (QD) 
is used to characterize the quantum correlations. 

Recently, many investigations have been made in un- 
derstanding the dynamical creation and evolution of cor- 
relations between the nearest-neighbor particle pairs and 
between two distant particles which are not connected by 
direct interactions in the spin chain model, for example, 
XX, XY , and Ising model systems @~Qj]|- As we know, 
if two particles directly interact with each other, then 
the correlations between these two particles can be built 
dynamically from an initial state without correlations. 



However, if two distant particles in a quantum network 
indirectly interact through local interactions, how fast 
will the correlations between these two distant particles 
be dynamically generated? The Lieb-Robinson theorem 
[12Hl7| gives an intriguing answer to this question: The 
speed of the correlations between two distant particles 
has an upper bound, which is called the Lieb-Robinson 
velocity. In other words, the correlations outside the 
light cone defined by the Lieb-Robinson velocity can be 
neglected. Recently, the Lieb-Robinson theorem has re- 
ceived renewed interest and has been applied to the con- 
densed matter theory and quantum information theory 
fl8U22j ]. For example, it can be used to derive a general 
relation on the two-party correlations in the many-body 
ground states. 

As an upper bound on the speed of correlation gener- 
ation, can the Lieb-Robinson velocity be observed in a 
concrete quantum network? As far as we know, there 
is no related experimental report so far. In fact, we 
can obtain the Lieb-Robinson velocity through finding 
the maximum generalized Lieb-Robinson velocity in all 
kinds of conditions. Here the generalized Lieb-Robinson 
velocity is defined as the correlation propagation speed 
from a given initial non-correlated state, which is upper 
bounded by the Lieb-Robinson velocity. In this paper, 
our central aim is to demonstrate how to observe the 
generalized Lieb-Robinson velocity in an Ising spin chain 
with a strong transverse field. To solve this problem, 
we first need to choose a model whose dynamics rela- 
tively easy to simulate. Then we need the measures to 
characterize different types of correlations. Finally, we 
need to give a criterion to judge whether the correla- 
tions between two distant particles appears. It should be 
pointed out that the generalized Lieb-Robinson velocity 
we study in this article is not the upper bound of speed 
but the concrete propagation velocity of correlation. In 
this article, we choose the transverse field Ising chain 
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as the basic model, which is exactly solvable for eigen 
problems. Wc consider the measures of different types of 
correlations, including correlation functions, mutual en- 
tropy EoF, and QD. Thus we can investigate whether the 
generalized Licb- Robinson velocity depends on the corre- 
lation measures. The criterion for correlation appearance 
is to set the correlation measure to a value numerically so 
small such that the correlation appearance time is almost 
fixed. 

In our article, we find that the generalized Lieb- 
Robinson velocity can be obtained by analyzing a corre- 
lation measure. Almost all types of correlation measures 
can demonstrate the generalized Licb-Robinson velocity 
in their dynamics. In particular, we find a sudden switch 
of the generalized Lieb-Robinson velocities with the in- 
creasing of the spin number. This paper is organized as 
follows. In Sec. [ill wc will briefly introduce the physical 
model and the solution of the dynamics for the reduced 
two-particle state. In Sec. IIII1 we will study and com- 
pare the dynamics of different kinds of correlations, and 
show the relation between the evolution of correlations 
and the length of the spin chain. In Sec. IIV| we will 
study the buildup of correlations. Finally wc will give a 
brief summary. 



II. MODEL, APPROXIMATION AND 
SOLUTION 

We will consider a transverse field Ising spin chain 
(TFIC) with free ends, whose Hamiltonian can be written 
as 

JV-l N 
t=l i=l 



A. Rotating-wave approximation for TFIC in a 
strong magnetic field 

In the weak coupling region (J -C B), the Hamiltonian 
(fljl can be rephrased in the following more convenient 
form 

w-i 

H = ~ J £ (°T +0 <+! + + °f + *<+"i 

N 

^"4i)- fl E<. ( 2 ) 

i=l 

where arc the two eigenvectors of erf, and 

°i = l/ u )i w ith fi,v € {+ . — } . Under the rotating- 
wave approximation (RWA) [23j , the Hamiltonian be- 
comes 

N-l N 

i=i i=i 

The Hamiltonian -ffrswA is unitarily equivalent to the 
XY spin model. Notice that in the above Hamiltonian 
([3]), the x-component of the total spin erf = J^i=i CT f i s 
conserved. In other words, an invariant subspace of the 
Hamiltonian can be characterized by a given eigenvalue 
of erf, and the dynamics of the system can be studied 
independently in these invariant subspaces. 

Before investigating further the dynamics of the sys- 
tem, we directly check the validity of RWA for the Hamil- 
tonian by numerically comparing the eigenvalues and the 
eigenvectors between the Hamiltonian without RWA (fT|) 
and the one with RWA Q. The numerical results for 
A = 3 are demonstrated in Fig. [2] As expected, wc 
find that RWA is an excellent approximation as long as 
S» J. 



where J is the coupling constant, B is the strength of 
the transverse filed, N is the total number of spin, and 
af (a = x, y, z) are the Pauli operators acting on the i-th 
spin of the chain. 

To investigate the dynamical creation of correlations 
along a TFIC, we consider the following physical pro- 
cess. First, let the system stay in the ground state of the 
system. Then flip the first spin. Our aim is to observe 
how the correlations between the first spin and the last 
spin will be created dynamically. The physical setting is 
demonstrated in Fig. [T] 




FIG. 1: Dynamical generation of correlations between the first 
spin and the last spin in the spin chain (TV = 7). 



B. Quantum state evolution under RWA 

In this section, we will obtain the analytical results on 
the quantum state evolution of the system under RWA. 

In the ground state of the Hamiltonian Hrwa, all 
the spins point along the positive direction of x axis. 
While flipping the first spin, we get the initial state 
|V(0)> = |- + + ---+). Because of#(0)) = (A-2)|V>(0)) 
and [erf ,_ffriWA] = 0, the quantum state will evolve with 
time in the eigen space of erf. with eigenvalue N — 2. All 
the eigenvectors of the subspace can be denoted as \m) = 
n™! 1 l+)i|-)™nf= m+ i \+h for m e {1,2, • • • ,N}. Us- 
ing these new notations, the initial state |^(0)) = |1), 
and the Hamiltonian in this subspace 
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FIG. 2: (Color online) (a) The energy spectrum (N — 3) as a 
function of B/J. Dashed lines: the result in RWA; Solid lines: 
the result without RWA. (b) The overlap between the wave 
function obtained with RWA and without RWA as a function 
of B/J. Dashed line: the overlap of ground states; Dotted- 
dashed line: the overlap of the first exctied states; Sloid line: 
the overlap of the second excited states. 



Since we will consider the correlations between the first 
spin and the last spin, we only need the reduced density 
matrix of these two spins, which is given by 



Pin = 









l^ll 2 

AiA* N 

\ o o 





A\A N 




(9) 



i-i^r-Mwiv 



The above formula implies that only A\ and An are 
needed to be calculated for our purpose, which reduces a 
large amount of computation. 



III. DYNAMICS OF DIFFERENT TYPES OF 
CORRELATIONS 

To study the correlations between the first spin and the 
last spin, we adopt both traditional method and informa- 
tion method to characterize the degrees of correlations. 
In this section, we will numerically study the dynamical 
evolution of different measures of correlations between 
the first spin and the last spin. The correlation measures 
we adopt include the correlation functions, the mutual 
information, quantum discord, and the entanglement of 
formation. We aim to find out the relation between the 
evolution of correlations and the length of spin chain. 



The Hamiltonian i/rwA has eigenstates 

n— 1 v ' 

with eigenvalues 

E k = -(N-2)B + 2Jcosj^, (6) 

for k = 1,2,- •• ,N. 

The quantum state of the system at time t is given by 

N 

|*(t)) = 5>„|n), (7) 

71=1 

where 



A. Traditional method: Correlation function 

Correlation function (CF) is a traditional tool in de- 
scribing the correlation effects in a many-body system. 
For a two-qubit system, the CFs are defined by [24| 

C F (o?,ofy = Tr (pia«) - Tr Tr (p N afy , 

(10) 

where a,/3 E {x,y, z}, and p\ and pjv are the reduced 
density matrices of the bipartite quantum state Pin- In- 
serting Eq. ([9]) into Eq. (fTU|). we have 

C F (ala z N ) = -ilArf \A N \ 2 , (11) 
CfK,4) = A\A* N + A\An- (12) 

Here we find that Cf (cf ,ct^-) = when N is even. 
This result can be proved as follows. We have the relation 
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C F {crl,a x N ) 



A 1 A* N + A* 1 A N 



N 



— £ cos[(E k -E m )t}(^l) m+1 sm 2 



(N+l) 



k,m— 1 
JV 



kn 
N+l 



sin 



7V + 1 



(7V+1) 2 



cos 



fc,m— 1 



2J cos 



iV + 1 
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kn 
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x(-l) 
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sin 



kn 
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N+l 
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cos 



;=1 fe=i 
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(13) 



o 

_ -0.2 

N 2 

D 

V-0.4 
°"-0.6 
-0.8, 



(a) 



^-0.01 
o u --0.02 

-0.03L 



n 



25 
t 



25 
t 



50 



50 



(d) 



V-0.005 
O 

-0.01„L 






( b) t| 








^_ .02 








r 










o u --0.04 








-o.oa 









25 
t 



50 




(b) 



0.4 
_ 0.2 

x Z 
D 

V 

LL 

° ^0.2 
-0.4, 



— 0.05 

x Z 
D 

V" 0- 

°-0.05 

-0.1 L 




(d) 



25 
t 



25 
t 



50 




50 



FIG. 3: Dynamics of correlation function Cf (erf, o^) of the 
spin chain with different lengths: (a) N = 10, (b) N = 20, 
(c) N = 30, (d) N = 40. 



FIG. 4: Dynamics of correlation function Cf (cf , <J%) °f t ne 
spin chain with different lengths: (a) N — 9, (b) N = 19, (c) 
N = 29, (d) N = 39. 



When TV is even, we have CV (of , cTjy ) = — Cp{a\ , c^-), 
which implies that Cf(v\ , crff) = 0. 

It is worth pointing out that Cp{a{^a z N ) 7^ whether 
iV is odd or even. These different correlation behaviors 
between Cp(af, c^y) and Cp(af, aff) shows that the cor- 
relations measured by correlation functions depend on 
which correlation function we choose. In the viewpoint of 
quantum information, correlation in a two-partite quan- 
tum state is a property of the state, and the characteri- 
zation of the correlation does not depend on what mea- 
surements we obtain the correlation. This information 
viewpoint will be detailedly discussed in the next section. 

The numerical results about the dynamical evolution 
of the two CFs, Cp (af,a^) and Cp (af,aff), of the spin 
chain with different lengths are shown in Fig. [3] and Fig. 
3J From the two figures, we can observe that the creation 



of correlations between the first spin and the last spin, 
which have no direct interactions, is not instantaneous 
but requires the formation time, and the formation time 
is directly proportional to the length of spin chain. In 
other words, the longer the distance between the first 
spin and the last spin is, the more time the appearance 
of their correlation needs. 



B. Information method: Mutual information, 
quantum discord, entanglement of formation 

In quantum information science, the characterization 
of correlations in a bipartite quantum state is relatively 
well understood. The degree of the total correlation in 
a bipartite quantum state is measured by the mutual 
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FIG. 5: Demonstration of classification of the total correlation 
meausred by mutual entropy. Classical correlation, quantum 
discord, quantum entanglement are abbraviated as CC, QD, 
and QE respectively. 



entropy. The total correlation can be classified into clas- 
sical correlation and quantum discord [1, [|[ . Quantum 
entanglement is a special type of quantum discord. The 
classification of correlations is demonstrated in Fig. \5\ 

The total correlation between two subsystems A and 
B arc quantified by the mutual information (MI) 



l{p AB ) = S{p A ) + S{p»)-S{p A ») 



AB\ 



(14) 



where S(p) = — Tr(plog 2 p) is the von Neumann entropy. 
To define the classical correlation contained in the state 
p( AB \ we consider the following process. B performs a 
projective measurement {Tlj} on the subsystem B, and 
will get a state 



1 



pt = — Tr B (7 ® U k p AB I ® n fe j 
Pk 



(15) 



with the corresponding probability p k = 
Ttab (7 ® Ii-kP AB I ® life) . By performing a measure- 
ment on subsystem A, A wants to give the information 
on which measurement result got by B. The upper 
bound of the information is the Holcvo bound 



For a two-qubit system, all the above correlation mea- 
sures can be obtained analytically or numerically. For 
example, in our case, the mutual information is 



1(1, N) = \A N \ 2 log 2 



\Ai\ 2 \og 2 



l-U 



N 



N\ 



l^ni-IArl 2 



1-4 



N 



(l-|Ar| 2 ) (\ Al \ 2 + \A N \ 



log 2 



\Ai\'[l 

l-l^l 2 -^ 



|Ai| 2 -L4 



A' 



AM 



l^ll 



1 - \A 



N 



(20) 



which is shown in Fig. 6. Following the method intro- 
duced in Ref. [25l [2q |. we can numerically obtain the 
dynamical evolution of QD as shown in Fig. 7. In or- 
der to compare the dynamical evolution between MI(QD) 
and EoF, we also plot the dynamical evolution of EoF as 
shown in Fig. 8. 
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X (p AB \{n k }) = S(p A )- Pk S(p A ). 



(16) 



Then a measure of classical correlation (CC) in the state 
p (AB ^ is defined by 1 



J{p 



AB \ 



max X (0{n fe }). 
{iiki 



(17) 



Once CC is obtained, the QD is obtained by subtracting 
CC from the MI 



V{p AB )=l{p A »)-J{p A "). 



AB\ 



„AB\ 



(18) 



Quantum entanglement is a special type of quantum 
discord, which is widely investigated in quantum infor- 
mation. One of the most useful measures of entanglement 
is entanglement of formation 0, 0] , which is defined by 

£{P [AB) ) = \ min Pi l{\4 AB) )), (19) 
where p^ = Eift I- 



FIG. 6: Dynamics of MI of the spin chain with different 
lengths: (a) N = 20, (b) N = 30, (c) N = 40, (d) N = 500. 

In Fig. [51 Fig. \7\ and Fig. [51 we plot four cases in 
each figure, that is, N = 20, N = 30, N = 40, and N = 
500, in order to find the relation between the dynamical 
evolution and the length of spin chain. From these three 
figures, we observe that the dynamical evolutions of MI, 
QD and EoF are similar with CFs, that is, the creation 
of MI, QD and EoF between the first spin and the last 
spin is not instantaneous but requires the formation time, 
and the longer the spin chain is, the more time to create 
correlations takes. In addition, the maximum amplitudes 
that the MI, QD and EoF can reach all become smaller 
with the increase of the length of spin chain. However, 
the degrees of amplitude reduction are not same. For 
example, we consider TV = 500, from Fig |6(d)| Fig |7(d)| 
and Fig |8(d)| we can see that the maximum amplitudes 
of MI and QD are about 4 x 10 -2 , while the maximal 
amplitude of EoF is only about 6 x 10 -3 , which implies 
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FIG. 7: Dynamics of QD of the spin chain with different 
lengths: (a) TV = 20, (b) TV = 30, (c) TV = 40, (d) TV = 500. 
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FIG. 8: Dynamics of EoF of the spin chain with different 
lengths: (a) TV = 20, (b) TV = 30, (c) TV = 40, (d) TV = 500. 



that MI and QD are more robust to the length increase 
of the spin chain than EoF. 



IV. START-UP TIME WITH DIFFERENT 
CHAIN LENGTHS 

In the above section we have shown that the time 
needed to create correlations between the first spin and 
the last spin increases with the increase of the spin chain 
length. However, we still do not know how fast correla- 
tions can be created between the first spin and the last 
spin in our concrete model, though there is an obvious 
and correct answer: No faster than the Lieb-Robinson 
speed. In this section, we will numerically study the re- 
lation between the start-up time of correlations for dif- 
ferent correlation measures and the length of spin chain. 
We aim to find out the propagation speed of correlations 
and investigate whether the propagation speed depends 



FIG. 9: (Color online) Start-up time of (a) MI, (b) EoF, (c) 
QD, and (d) CC as a function of TV with different criteria. 
The criteria of solid line, dashed line, and dotted-dashed line 
are 8 = 10~ 4 , 8 = 10 -5 , and 8 — 10 -6 , respectively. 



on the correlation measures. 

Fig. [9] shows the start-up time of MI, QD, EoF, and 
CC as a function of TV, the length of spin chain, where TV 
is an integer we choose from 2 to 1000. Here we define the 
criterion for correlation appearance, which can be used 
to judge whether the correlations have been created or 
not at a given time. In every subfigure of Fig. [HJ we 
plot three curves. Different curves have different criteria 
of start-up 5, that is, the criteria of solid line, dashed 
line, and dotted-dashed line are S = 10~ 4 , S = 10~ 5 , and 
S = 10~ 6 , respectively. From Fig. [9j we observe that the 
start-up time as a function of the spin chain length has 
two segments, and there is a sudden switch point, which 
moves rightwards with the decrease of the criterion. From 
Fig. |6(d)[ Fig. |7(d)| ,and Fig. |8(d)[ we can see that there 
are a series of obvious peaks in the figures. In fact, the 
first segment in Fig. [9]appears since the first peak reaches 
the criterion, and the second segment appears since the 
second peak reaches the criterion. Thus we plot the first 
peak value of MI as a function of the spin chain length as 
shown in Fig. 10(a) We find the relation between log I 
and log TV is linear, that is 



I = N a e , 



(21) 



where a = -2.8529 and (3 = 3.0981. We also plot the 
second peak value of MI as a function of the spin chain 
10(b)| the relation between logX 
3 ut a = -0.9257 and j3 = 0.5725. 



length as shown in Fig. 
and log TV is also linear, 



So the first peak will disappear sooner than the second 
peak with the increase of TV. This implies that the first 
peak becomes lower than the criterion with the increase 
of TV, but the second peak is still higher than the cri- 
terion. This is the reason for the appearance of switch 
point in Fig. [51 

Lastly, in order to investigate whether the propagation 
speed depends on the correlation measures, we plot them 
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FIG. 10: log I is linearly proportional to log N for (a) the first 
peak value, (b) the second peak value. 



in the same figure with the criterion S = 10 -6 as shown 
in Fig. [11] to compare the start-up times of MI, QD, 
CC and EoF. We find that the propagation velocity does 
not depend on the correlation measures, but the sudden 
switch point depends on the correlation measures. 




exponential indexes. It is these different characteristics 
of correlations that lead to the sudden switch of the Lieb- 
Robinson velocities. 



DISCUSSION AND CONCLUSION 



We have numerically investigated the problem on the 
generalized Lieb-Robinson velocities in TFIC. We show 
that the generalized Lieb-Robinson velocity can be ob- 
served by the dynamics of the correlations between the 
first spin and the last spin. A sudden switch of the gen- 
eralized Lieb-Robinson velocities appears in TFIC. Here 
a main open problem is whether the phenomena in TFIC 
are model independent. However, the method we adopt 
for TFIC uses the analytical results, which is not avail- 
able for a general spin model. In addition, as far as we 
know, the algorithms for quantum dynamics are not pow- 
erful enough to simulating for a sufficient long time with 
a long spin chain. In this direction, we study another 
related mode, the isotropic Heisenbcrg chain. 

The start-up time of CF of isotropic Heisenbcrg chain 
with uniform coupling as a function of the chain length N 
is demonstrated in Fig. [TSJ We find that the phenomena 
on the generalized Lieb-Robinson velocity have the same 
characteristic in the isotropic Heisenberg chain. 
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FIG. 12: Start-up time of CF of the Heisenberg chain with 
uniform coupling as a function of N: (a) from N = 2 to 
N = 200, (b) from N = 614 to N = 637. 



FIG. 11: (Color online) Comparison of start-up times of MI, 
QD, CC and EoF with the criterion 5 = 10" 6 . 

In the above numerical results, the correlation start- 
up time is linearly proportional to the length of the 
spin chain, which implies the existence of the general- 
ized Lieb-Robinson velocity. The sudden switch of the 
generalized Lieb-Robinson velocities implies that differ- 
ent types of correlations may have different generalized 
Lieb-Robinson velocities. In the present model, the am- 
plitudes of these two types of correlations decreases ex- 
ponentially with the length of the spin chain, and these 
two types of correlations are characterized by different 



In conclusion, we have numerically investigated the 
generalized Lieb-Robinson velocity in an Ising spin chain 
with a strong transverse field through studying the dy- 
namical evolution of correlations. The generalized Lieb- 
Robinson velocities are demonstrated in different types of 
correlation measures, which include correlation function, 
mutual information, quantum discord, and entanglement 
of formation. We find that one of the correlation func- 
tions shows a special behavior depending on the parity 
of the spin number. In particular, we find that there is a 
switch in the generalized Leb-Robinson velocities, which 
implies that different types of correlations may show dif- 
ferent generalized Lieb-Robinson velocities. 
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